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NASA-CR-202613 ABSTRACT 

The origin of the highly eccentric, inclined, and resonance-locked orbit of Pluto has long been a puzzle. A 
possible explanation has been proposed recently [Malhotra, 1993, Nature, 365, 819] which suggests that 
these extraordinary orbital properties may be a natural consequence of the formation and early dynamical 
evolution of the outer solar system. A resonance capture mechanism is possible during the clearing of the 
residual planetesimal debris and the formation of the Oort Cloud of comets by planetesimal mass loss from 
the vicinity of the giant planets. If this mechanism were in operation during the early history of the planetary 
system, the entire region between the orbit of Neptune and approximately 50 AU would have been swept 
by first-order mean motion resonances. Thus, resonance capture could occur not only for Pluto, but quite 
generally for other trans-Neptunian small bodies. Some consequences of this evolution for the present-day 
dynamical structure of the trans-Neptunian region are (i) most of the objects in the region beyond Neptune 
and up to —50 AU exist in very narrow zones located at orbital resonances with Neptune (particularly the 
3:2 and the 2:1 resonances); and (ii) these resonant objects would have significantly large eccentricities. The 
distribution of objects in the Kuiper Belt as predicted by this theory is presented here. 


1. INTRODUCTION 

In the widely accepted paradigm for the formation of the 
solar system, the planets accumulated in a highly dissipative 
disk of dust and gas orbiting the protosun, and most planets 
formed in near-circular and nearly coplanar orbits. The out- 
ermost planet, Pluto, is an oddity in this scheme; its orbit is 
highly eccentric (e = 0.25) and inclined (17° to the ecliptic). 
The large eccentricity means that Pluto crosses the orbit of 
Neptune, and it traverses a very large region of space from 
just inside the orbit of Neptune at 30 AU to almost 50 AU 
Indeed, soon after Pluto was first discovered in 1930, it was 
realized that the dynamical lifetime of this new planet was 
short before a close encounter with Neptune radically altered 
its orbit (Lyttleton 1936). Three decades later, it was found 
that a dynamical protection mechanism exists that prevents 
close encounters between Pluto and Neptune: a 120,000 yr 
orbit integration of the outer planets by Cohen & Hubbard 
(1965) showed that Pluto is locked in a 3:2 orbital resonance 
with Neptune which maintains a large longitude separation 
between the planets at orbit crossing and causes Pluto’s peri- 
helion to librate about a center ±90° away from Neptune. 
Since the Cohen and Hubbard work, orbit integrations of 
increasingly longer times have uncovered several other reso- 
nances and near-resonances in Pluto’s motion (Williams & 
Benson 1971; Applegate et al. 1986; Sussman & Wisdom 
1988; Milani et al. 1989). Perhaps the most important of 
these “weaker” resonances is the “argument-of-perihelion 
libration” which ensures that at perihelion Pluto is close to 
its maximum excursion above the mean plane of the solar 
system; this has the effect of increasing the minimum ap- 
proach distance between Pluto and Neptune and between 
Pluto and Uranus than would otherwise be the case. (It is 
worth emphasizing that the longitude-of-perihelion libration 


centered at ±90° relative to Neptune and the argument-of- 
perihelion libration about 90° are two quite distinct phenom- 
ena; the former is associated directly with the 3:2 orbital 
period resonance lock between Neptune and Pluto.) Pluto’s 
orbit is confined in a very narrow region of relative orbital 
stability near the 3:2 Neptune resonance, a region in phase 
space that is bounded by highly chaotic orbits. (See Malhotra 
& Williams 1994 for a recent review of Pluto’s orbital dy- 
namics). 

There has, of course, been much speculation as to the 
origin of this extraordinary orbit of Pluto (Lyttleton 1936; 
Harrington & Flandem 1979; Farinella et al. 1979; Dormand 
& Woolfson 1980; Olsson-Steel 1988; Malhotra 1993a; Le- 
vison & Stem 1995), all but one — Malhotra (1993a, hereaf- 
ter referred to as Paper I) — of these speculations requiring 
one or more low-probability “catastrophic” events. A popu- 
lar theme in the earlier scenarios was that Pluto is an escaped 
satellite of Neptune. On the other hand, studies and modeling 
of the physical characteristics and composition of Pluto have 
led to the conclusion that it accumulated in a heliocentric 
orbit of ice-rich planetesimals in the outer reaches of the 
Solar Nebula rather than in a circumplanetary disk (McKin- 
non & Mueller 1988; Tancredi & Fernandez 1991). In con- 
sonance with this, the recent theories by Malhotra (1993a) 
and Levison & Stem (1995) both propose that Pluto formed 
in a near-circular coplanar heliocentric orbit beyond the or- 
bits of the giant planets, but they differ in the physical and 
dynamical mechanisms that placed Pluto in its unusual orbit. 

In the model proposed by Levison and Stem, the orbital 
configuration of the planets is taken as observed today, ex- 
cept that a “test Pluto” is placed in an initially low- 
eccentricity, low-inclination orbit near the 3:2 Neptune reso- 
nance. With some fine tuning of initial conditions, orbit 
integrations show that such an orbit has its eccentricity and 
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inclination pumped up to values comparable to those of the 
real Pluto in a time scale of about 10 7 yr. However, the orbit 
remains chaotic during this evolution; Levison and Stern 
then propose that Pluto was “knocked” into the stable 3:2 
resonance libration region by one or more dissipative colli- 
sions with a neighboring small body or bodies. 

Malhotra’s theory (Paper I) does not invoke any cata- 
strophic collisions, and is possibly compatible with the stan- 
dard paradigm for planet formation. [The reader is referred to 
Levy & Lunine (1993) for reviews of planet formation 
theory.] In this model, an initially low-inclination, nearly cir- 
cular orbit of Pluto beyond the orbits of the giant planets 
evolves into its Neptune-crossing but resonance-protected or- 
bit as a result of early dynamical evolution of the outer solar 
system. The physical causes of this evolution lie in the late 
stages of planet formation when the gravitational scattering 
and eventual clearing of remnant planetesimal debris by the 
giant planets (and the concomitant exchange of energy and 
angular momentum between the planets and the planetesi- 
mals) may have caused a significant evolution of the giant 
planet orbits (Fernandez & Ip 1984). In particular, Neptune’s 
orbit may have expanded considerably, and its exterior or- 
bital resonances would have swept through a large region of 
trans-Neptunian space. During this resonance sweeping, 
Pluto could have been captured into the 3:2 orbital period 
resonance with Neptune and its eccentricity (as well as 
inclination — see Sec. 4) would have been pumped up during 
the subsequent evolution. 

The dynamical mechanisms invoked in this theory are 
quite general, and would apply not only to the evolution of 
the trans-Neptunian body labeled “Pluto,” but also to any 
other members of the trans-Neptune region. While this pos- 
sibility was implicit in Paper I, it is my purpose in the 
present paper to make explicit the implications and predic- 
tions of this “resonance capture theory” of the origin of Plu- 
to’s orbit for the present-day architecture of the Solar system 
beyond Neptune. 

That the outermost parts of the Solar system may be 
populated by primordial icy planetesimals has been conjec- 
tured on both theoretical and observational grounds. For ex- 
ample, Kuiper (1951) suggested this on the basis of theoreti- 
cal considerations of the genesis of the planetary system 
from the primordial Solar Nebula. Whipple (1964) and 
Bailey (1983) speculated on a massive comet belt as the 
source of unexplained perturbations of Neptune’s orbit [al- 
though this argument must now be discarded as the post- 
Voyager revisions in the planetary ephemeris no longer show 
any unexplained residuals in Neptune’s motion (Standish 
1993)]. Hamid et al (1968) analyzed the orbital plane per- 
turbations of comet P/Halley and concluded that any comet 
belt between 40 and 50 AU has a total mass less than L-#® . 
More recently, it has been suggested that the observed short- 
period comets with orbital periods ^20 yr originate in a belt 
of low- inclination bodies just beyond the orbit of Neptune, 
between 35 and 50 AU (Fernandez 1980; Fernandez & Ip 
1983). The older hypothesis that short-period comets origi- 
nate in a population of near-parabolic Oort Cloud comets 
(which are perturbed into shorter orbits by the giant planets) 
appears unlikely: Duncan et al (1988), Quinn et al (1990) 


have shown that the orbital element distribution of the ob- 
served short-period comets is inconsistent with a source in 
the nearly isotropic Oort Cloud but is compatible with a 
disk-like source in a trans-Neptune comet belt, which they 
call the “Kuiper Belt.” A possible member of the Kuiper 
Belt was first discovered in 1992 at a distance of 41 AU from 
the Sun (1992 QBj, reported in Jewitt & Luu 1993), and 
several additional discoveries have been reported since (Je- 
witt & Luu 1995). 

The dynamical structure of this putative comet population 
as determined by the long-term (conservative) gravitational 
perturbations by the planets has been the subject of two re- 
cent studies (Levison & Duncan 1993; Holman & Wisdom 
1993). These studies assumed a uniformly distributed initial 
population in near-circular, low-inclination orbits, and the 
planets in their present orbital configuration; they sougnt to 
determine the extent and nature of any orbital instabilities 
that might operate on billion year time scales to transport 
those putative comets into planet-crossing orbits. However, 
any trans-Neptune population of planetesimals was undoubt- 
edly subject to dynamical evolution during the planet forma- 
tion process, and the initial conditions assumed in the above 
studies are not necessarily representative of the state of the 
Kuiper Belt at the end of planet formation, as acknowledged 
in Holman & Wisdom (1993). Here I discuss the “dynamical 
sculpting” of the Kuiper Belt that would have occurred due 
to the early orbital evolution of the outer planets (during the 
late stages of their formation) as predicted by the “resonance 
capture theory” for the origin of Pluto’s orbit. The results of 
this study indicate that the Kuiper Belt would have been 
“sculpted” into a highly nonuniform distribution early in 
solar system history, and this structure would be largely pre- 
served to the present epoch: the region beyond Neptune’s 
orbit and up to approximately 50 AU heliocentric distance 
should have most of the primordial small bodies locked in 
orbital resonances with Neptune, particularly the 3:2 and the 
2:1 orbital resonances which are located at semimajor axes 
of approximately 39.4 and 47.8 AU, respectively. 

The rest of this paper is organized as follows. In Sec 2, I 
summarize the “resonance capture theory” for the origin of 
Pluto’s orbit. In Secs. 3 and 4, 1 describe the numerical simu- 
lations conducted to determine the implications of this theory 
for the dynamical structure of the trans-Neptunian solar sys- 
tem. In Sec. 5, I provide a discussion of the results, a com- 
parison with previous theoretical studies of the Kuiper Belt, 
and some consequences for observational surveys of the 
outer solar system. 

2. RESONANCE CAPTURE THEORY FOR THE ORIGIN 
OF PLUTO’S ORBIT 

The phenomenon of capture into resonance as a result of 
some slow dissipative forces is common in nature, and there 
exists a large body of literature devoted to its study. A well- 
developed solar system example is the formation of orbit- 
orbit resonances among the satellites of the giant planets by 
the action of slow tidal dissipation (see Peale 1986 and Mal- 
hotra 1994a for reviews). In general, capture into a stable 
(long-lived) orbit-orbit resonance is possible when the orbits 
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of two bodies approach each other as a result of the action of 
some dissipative process. In Paper I, I proposed that Pluto 
may have been captured into the 3:2 resonance with Neptune 
during the late stages of planet formation, when Neptune’s 
orbit expanded outward as a result of angular momentum 
exchange with residual planetesimal debris. The physics of 
this mechanism is summarized below. 

The giant planets’ gravitational perturbations were instru- 
mental in clearing their interplanetary regions of the residual 
unaccreted planetesimal debris. It is believed that the Oort 
Cloud — which is a roughly isotropic distribution of comets 
surrounding the planetary system at distances in excess of 
~10 4 AU — was populated by icy planetesimals scattered 
outward from the vicinity of the giant planets; estimates of 
the total mass of the Oort Cloud are in the range 10- 10 2 
, based upon observations of long period comets and 
extensive theoretical modeling (cf. Weissman 1990). While 
the formation and dynamical evolution of the Oort Cloud has 
been a subject of extensive research (indeed, it is subspe- 
cialty within planetary science), relatively little attention has 
been given to the back reaction on the planets themselves of 
the planetesimal scattering process that populated the Oort 
Cloud. 

Consider the scattering of a planetesimal of mass m c by a 
planet of mass at orbital radius a . If the planetesimal is 
initially in a near-circular orbit similar to that of the planet 
and is ejected to a solar system escape orbit, it follows from 
conservation of angular momentum that the planet suffers a 
loss of orbital angular momentum and a corresponding 
change of orbital radius, Sa , given by 1 


For a planetesimal scattered outward, but one that does not 
achieve a solar system escape orbit (remaining bound, for 
example, in the Oort Cloud), the numerical coefficient on the 
right-hand side of Eq. (1) would be slightly smaller than 
unity. Conversely, planetesimals scattered inward would 
cause an increase of orbital radius and angular momentum of 
the planet. A single, massive planet scattering a population of 
planetesimals in near-circular orbits in the vicinity of its own 
orbit would, to first order, suffer no net change of orbital 
radius as it scatters approximately equal numbers of plan- 
etesimals inward and outward. 

However, the four Jovian planets acting together evolve 
differently from this simple picture, as first pointed out by 
Fernandez & Ip (1984) who modeled the late stages of ac- 
cretion of planetesimals (“proto-comets”) by the proto-giant 
planets (and the concomitant exchange of energy and angular 
momentum between the planetesimals and the planets). Their 
numerical simulations showed a small decrease in orbital 
radius for Jupiter and significant increases in orbital radius 
for Saturn, Uranus, and Neptune. The reason for this orbital 
evolution — in particular that of Neptune — may be under- 
stood by means of the following heuristic picture of the 

lr The expression given in Eq. (1) in Malhotra (1993a) is incorrect in that the 
coefficient on the right-hand side of the equation was underestimated by 
about 20%; I am indebted to S. Tremaine for calling my attention to this. 
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Fig. 1. A schematic diagram to illustrate the outward radial migration of 
Neptune and its exterior orbital resonances during the late stages of planet 
formation. The distance from the Sun is along the horizontal direction. Nep- 
tune’s outward orbital migration is shown along the path marked N-N. For 
clarity, only two first-order resonances (3:2 and 2:1) are shown (dotted 
lines). A “Pluto” in an initially circular, nonresonant orbit beyond Neptune 
could have been captured into the 3:2 resonance and would evolve along the 
solid line path indicated by P — 3/2. 

clearing of a planetesimal swarm from the vicinity of Nep- 
tune. Suppose that the mean specific angular momentum of 
the swarm is initially equal to that of Neptune. At first, a 
small fraction of planetesimals is accreted as a result of 
physical collisions, and of the remaining, there are approxi- 
mately equal numbers of inward and outward scatterings. To 
first order, these cause no net change in Neptune’s orbit. 
However, the subsequent fate of the inward and outward 
scattered planetesimals is not symmetrical. Most of the in- 
wardly scattered objects enter the zones of influence of the 
inner Jovian planets (Uranus, Saturn, and Jupiter). Of those 
objects scattered outward by Neptune, some are lifted into 
wide, Oort Cloud orbits while others return to be accreted or 
rescattered; a fraction of the latter is again (re)scattered in- 
ward where the inner Jovian planets control the dynamics. In 
particular the massive Jupiter is very effective in causing a 
systematic loss of planetesimal mass by ejection into solar 
system escape orbits. Therefore, as Jupiter preferentially re- 
moves the inward scattered Neptune planetesimals, the plan- 
etesimal population encountering Neptune at later times is 
increasingly biased toward objects with specific angular mo- 
mentum (and energy) larger than Neptune’s. Encounters with 
this planetesimal population produce effectively a negative 
drag on Neptune which results in Neptune experiencing a net 
gain of angular momentum and energy, hence an increase in 
its orbital radius. Note that Jupiter is, in effect, the source of 
this angular momentum and energy; however, owing to its 
much larger mass, its orbital radius decreases by only a small 
amount. 

If the above phenomenon did occur in the late stages of 
planet formation, one consequence of Neptune’s orbital ex- 
pansion is that its orbital resonances would have swept 
across a range of heliocentric distances comparable to its 
radial migration (see Fig. 1). During this “resonance sweep- 
ing,” a small body such as Pluto, initially in a near-circular 
orbit beyond Neptune, could have been captured and locked 
into an orbital resonance. 

Resonance capture is a complicated dynamical process 
and a very active subject of research in the nonlinear dynam- 
ics literature (see, for example, Wiggins 1991). The transi- 
tion from a nonresonant to a resonant orbit depends sensi- 
tively upon initial conditions, and in the context of the solar 
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system, upon the nature of the resonant as well as nonreso- 
nant gravitational perturbations, and the rate of orbit evolu- 
tion due to the dissipative effects. Under certain idealized 
conditions (“single resonance”), and in the limit of slow 
“adiabatic” orbit evolution, the probability of resonance cap- 
ture is relatively straightforward to calculate (Henrard & Le- 
maitre 1983; Borderies & Goldreich 1984). Such a calcula- 
tion shows that the capture probability for the 3:2 Neptune 
resonance is 100% for initial eccentricity (before the reso- 
nance encounter) less than —0.03; the capture probability 
decreases monotonically for higher initial eccentricities: it is 
less than 10% for initial eccentricities exceeding 0.15. 

Once an object was captured into an orbital resonance, 
perturbations from Neptune would transfer sufficient angular 
momentum to it to maintain it in the resonance by expanding 
its orbit in concert with that of Neptune. A byproduct of this 
evolution would have been the rapid excitation of the ob- 
ject’s orbital eccentricity. This is most readily seen in the 
following simplified analysis of the first-order perturbations 
of Neptune on a test-particle orbit. Close to an exterior 
/+!:; orbital resonance, the first-order perturbation equa- 
tions for the mean motion, n, and eccentricity, e, of the par- 
ticle are (cf. Brouwer & Clemence 1961) 

n = 3(j+ l)fi^n 2 ef(a)sin <fi y 

( 2 ) 

e = - sin <£, 

where jjb^=m^LS G is the mass of Neptune relative to the 
Sun, a= a n / a < 1 is the ratio of semimajor axes of Nep- 
tune and the test particle, f(a) is a positive function that can 
be expressed in terms of Laplace coefficients. 
<f>=(j+ 1)\— ;\ N -nr is the critical resonance angle, with X 
and xn the mean longitude and the longitude of perihelion, 
respectively, of the test particle, and \ N the mean longitude 
of Neptune. 

If the test particle is captured into resonance, its mean 
motion becomes locked to that of Neptune, so that the fol- 
lowing conditions hold: 

(3) 

(j+\)(n)=*j(n N ), 

where(n N ) is the rate of change of Neptune’s mean motion 
as its orbit expands. It then follows from Eqs. (2) and (3) that 

«n\ _ _J _/ 

*N/ 0'+l)\«N 

where the last equality follows from the Keplerian relation 
between mean motion and semimajor axis (/i 2 a 3 =const). 
Therefore, upon capture into resonance, the test particle’s 
eccentricity is pumped up at a rate determined by the average 
rate of expansion of Neptune’s orbit. The previous equation 
can be integrated to yield 


de 2 

~dt 


3(7 + 1) 



In 


a N .final 
fl N, initial 


(5) 


where refers to the value of Neptune’s semimajor 

axis at the point of resonance capture. Note that for initially 
near-circular orbits, the final eccentricity is insensitive to the 


initial eccentricity and depends only upon the extent of or- 
bital migration of Neptune. 

Applying this result to Pluto, one concludes that if an 
initially near-circular Pluto was captured into the 3:2 reso- 
nance 0 = 2) with Neptune and its eccentricity was reso- 
nantly pumped up to its current value of 0.25 in the subse- 
quent evolution, then Neptune was at a N «=* 25 AU at the 
point of resonance capture. Pluto’s initial orbital radius can 
then be inferred to have been near 33 AU The 5 AU outward 
radial migration inferred for Neptune must be regarded as a 
lower limit, as a N — 25 AU is the stage at which Pluto is 
inferred to have been captured in the 3:2 Neptune resonance. 
Other trans-Neptune objects captured in the 3:2 resonance at 
earlier or later times would have their eccentricities pumped 
up to larger or smaller values, respectively, than that of 
Pluto. Other first-order resonances of importance in trapping 
trans-Neptune objects are the 4:3 and the 2:1, currently lo- 
cated at semimajor axes of 36.5 and 47.9 AU, respectively. 
The 5:3 second-order resonance located at 42.3 AU also has 
a significant capture probability, and capture is also possible 
in other higher-order resonances such as the 7:5, and 7:4 
located at 37.7 and 43.7 AU, respectively. 

3. NUMERICAL MODEL FOR THE EVOLUTION 
OF THE TRANS-NEPTUNE POPULATION 

In the late stages of planetary accumulation, the exact 
magnitude of the radial migration of the Jovian planets due 
to their interactions with residual planetesimals is difficult to 
estimate without a full-scale A-body model. The work of 
Fernandez & Ip (1984) is suggestive, but simulations of this 
process stand to profit by the refinements in computal tonal 
technology that have occurred in the last ten years. Fcr ex- 
ample, the limitations of the software and hardware ava lable 
to Fernandez & Ip limited the total number of bodies in their 
simulations to about 2000; therefore, in order to start with a 
reasonable total initial mass in the planetesimal disk, the 
masses of their individual planetesimals were in the lather 
exaggerated range of (0.02-0. . Another approxima- 
tion in their model is the neglect of all but very close en- 
counters between the (proto-) giant planets and the planetesi- 
mals. I expect to improve upon this in future work. 

For the present, I do not attempt such detailed moceling 
here, but rather continue to use the simple model outlined in 
Paper I. Accordingly, the system consists of the Sun and the 
four Jovian planets with their present masses, together with a 
population of massless “test particles” representing rans- 
Neptunian objects. The massive planets perceive their full 
mutual gravitational interactions and also perturb the test 
particles, but the latter are noninteracting and do not perturb 
the planets. (It is perhaps worth emphasizing that there is 
little overlap between the trans -Neptunian population of 
small objects considered here and the population of residual 
planetesimals in the immediate vicinity of the giant planets 
whose gravitational interactions with the planets are pre- 
sumed to be driving the radial migration of those planets.) 
The orbital migration of the Jovian planets is modeled by a 
time variation of their semimajor axes according to th^ fol- 
lowing prescription: 
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a(t) = a f -ka exp(-f/r), (6) 

where df is the semimajor axis at the current epoch and 
Aa is the semimajor axis at the starting point (f — 0) 
of the simulation. (The epoch “/ = 0” refers to a time in the 
late stages of the genesis of the Solar system when the for- 
mation of the gas giant planets was largely complete, the 
Solar Nebula had lost its gaseous component, and the subse- 
quent evolution was dominated by gravitational interactions 
among the planets and the residual planetesimals. See Levy 
& Lunine (1993) for details on the various stages of solar 
system formation.) The orbit evolution given by Eq. (6) was 
implemented in the equations of motion by means of an ad- 
ditional “drag” force on each planet along the direction, v, 
of the orbital velocity given by 


Ar= 


v 

T 




(7) 


The numerical method described in Malhotra (1994b) was 
used for the orbit integrations; this method is based upon a 
second-order symplectic map (see Wisdom & Holman 1991) 
but modified for additional nongravitational forces. Note that 
the integration follows the orbits of the four major planets 
self-consistently (i.e., their mutual gravitational interactions 
are fully accounted for even as their orbits expand). 


4. RESULTS 

Here I report the results of several runs based upon the 
model described above. In all cases, the parameters for the 
planetary orbit evolution were as follows. The A a [cf. Eq. 
(6)] were chosen to be— 0.2, 0.8, 3.0, and 7.0 AU for Jupiter, 
Saturn, Uranus, and Neptune, so that the initial semimajor 
axes of these planets were approximately 5.4, 8.7, 16.3, and 
23.2 AU, respectively. These are within the ranges of radial 
displacements of these planets found in the Fernandez & Ip 
(1984) calculations. The planet masses and other initial or- 
bital elements were taken from Nobili et al (1989). 

The first two runs described here were aimed at determin- 
ing the current state of a primordial population of small ob- 
jects beyond Neptune in the Kuiper Belt (up to approxi- 
mately 50 AU heliocentric distance). In each of these runs, 
there were 120 test particles — representing the Kuiper Belt 
objects — with initial semimajor axes distributed uniformly in 
the range 28-52 AU and all angles (longitude of perihelion, 
longitude of ascending node, and mean longitude) chosen 
randomly from the range (0,27 t). In Run 1 (a “thin disk”), 
the initial eccentricities and inclinations were set to 0.01; in 
Run 2 (a “thicker disk”), the initial eccentricities and incli- 
nations were set to 0.05. The time scale r for the radial 
migration of the planets was taken to be 2X10 6 yr. The sys- 
tem was integrated for a period of 2X10 7 yr which is ten 
times the assumed orbital evolution time scale, r. At the end 
of the integration, the final orbits of the planets are very 
similar to their presently observed orbits. The integration of 
those test particles that suffered close approaches, i.e., within 
a Hill sphere radius, with any planet was terminated: the 
object was presumed to have failed to survive as a Kuiper 
Belt object. It was found that all objects with initial semima- 
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Fig. 2. The initial and final orbits of the surviving Kuiper Belt objects in 
Run 1 (“thin disk”) initial e = i = 0.01 shows by open circles; final ele- 
ments shown by triangles: and Run 2 (“thick disk”) initial e = / = 0.05 
show by solid circles; final elements shown by squares). For clarity, I show 
the final semimajor axes averaged over the last two million year period in 
the runs. The locations of several major orbital resonances with Neptune are 
indicated at the top of the upper panel. The concentrations of objects at these 
resonances have significantly large eccentricities and inclinations. 


jor axes less than 30 AU and a small fraction of those with 
initial semimajor axes between 30 and 34 AU failed to sur- 
vive. The survival rate in the “thicker disk” population was 
slightly lower than in the “thin disk” (85% vs 91%). (In 
general, the survival rate can be expected to be sensitive also 
to the planetary orbital evolution time scale, r, as discussed 
below.) Figures 2 and 3 summarize the results obtained in 
these two runs. Noteworthy features of the surviving test 
particle population are (i) the final semimajor axes are larger 
than 36 AU; (ii) the population is highly concentrated at 
primarily two resonances with Neptune — the 3:2 and the 2:1 
resonances located near 39.4 and 47.8 AU, respectively; 
smaller concentrations are also found at the 5:3, 4:3, 7:5, and 
7:4 resonances located at 42.3, 36.5, 37.7, and 43.7 AU, 
respectively; and (iii) the objects surviving in resonances 
have significant orbital eccentricities, typically 0.1 -0.3 (thus 
their perihelion distances can be as small as ~27 AU). 

Three other runs were made with a special focus on Pluto- 
like orbits captured in the 3:2 orbital resonance with Nep- 
tune. In each of these, 120 test particles with initially near- 
circular, coplanar orbits (e — i = 0.01) distributed uniformly 
in a e (29,35) AU were integrated for a period of 100 mil- 
lion years with the same model described above. The three 
runs differed only in the time scale rof orbit evolution of the 
planets; r was set equal to 2, 4 and 10 million years in Run 
3, Run 4, and Run 5, respectively. The survival rates in these 
runs were 79%, 68%, and 53%, respectively. As expected by 
design, most of the surviving bodies were locked in the 3:2 
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Fig. 3. The distribution of orbital semimajor axes, eccentricities, and incli- 
nations of the surviving Kuiper Belt objects in Runs 1 (lighter shading) and 
2 (darker shading). A large fraction of the objects beyond Neptune and up to 
a™* 50 AU are found in narrow resonance regions, as indicated in the top 
panel. (The distributions for Runs 1 and 2 are quite similar). 


resonance at the end of the evolution. (The exceptions were a 
few objects which were found locked in other nearby reso- 
nances.) An example of the typical orbital evolution of a test 
particle captured in the 3:2 resonance is shown in Fig. 4 
which displays the time variation of those variables that are 
of particular interest in the dynamics of Pluto’s orbit. Ob- 
serve that the semimajor axis stabilizes at the 3:2 resonance 
value; the eccentricity and inclination are both amplified; the 
3:2 resonance angle, 3\— 2\ N — uj (where \ N and X are 

the mean longitudes of Neptune and the test particle, and vj 
is the longitude of perihelion of the test particle), settles into 
stable libration about 180°, and the argument of perihelion, 
a>, also exhibits librations. Long-lived librations of o> about 
either +90° or -90° (i.e., perihelion near the maximum ex- 
cursion above or below the mean plane of the solar system) 
are found in approximately one-third of the objects trapped 
in the 3:2 Neptune resonance; often, the o) libration center 
hops between +90° and —90° over 10 7 yr time scales (cf. 
Fig. 4). The co librations are strongly correlated with large 
perturbations of the inclination. (This suggests that the reso- 
nance represented by the argument-of-perihelion libration 
may be responsible for the inclination excitation. However, 
this conjecture must await further analysis.) 

The distributions of the orbital elements of the survivors 
in these three runs are shown in Figs. 5-7. Noteworthy fea- 
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Fig. 4. A typical example of the orbital evolution of a Pluto-like Kuiper Belt 
object captured in the 3:2 resonance with Neptune. The 3:2 resonance angle, 
0=3\-2 X n -tb, settles into stable librations about 180°. There ensues a 
secular transfer of angular momentum from Neptune to the object that main- 
tains it in the resonance; the object’s orbital semimajor axis increases (in 
concert with Neptune’s) and its eccentricity is pumped up [cf. Eq. (2)]. In 
many cases, the inclination is also pumped up, but this depends sensitively 
upon the initial conditions of the orbit; the inclination behavior is highly 
correlated with the behavior of the argument of perihelion, io , which often 
exhibits long periods of libration about either +90° or -90°. 


tures are as follows. The eccentricity distribution is virtually 
identical in the three runs, with the vast majority of the ob- 
jects having e in the range 0.2 to 0.3 (as expected by design). 
Most of the objects remain in relatively low inclination 
(<10°) orbits, but a small fraction (up to 10% in Run 7) have 
their inclinations pumped up to higher values (15-20°), 
comparable to that of Pluto. There is a correlation between 
higher inclinations and larger r(i.e., slower evolution rate for 
the planetary orbits). 

5. CONCLUSIONS AND DISCUSSION 

The unusual properties of Pluto’s orbit may be a natural 
consequence — and a signature — of the early dynamical evo- 
lution in the outer solar system. The studies presented here, 
together with those reported in previous work Malhotra 
(1993, 1995) provide support for this case. I had shown pre- 
viously that if Neptune’s orbit expanded during the late 
stages of planet formation, then Pluto could have beer cap- 
tured from an initially near-circular, low inclination norireso- 
nant orbit beyond Neptune, into its current 3:2 resonant or- 
bit; during this evolution, its eccentricity would have been 
pumped up to the observed Neptune-crossing value while 
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Fig. 5. The final distribution of a, e, and / in Run 3 in which the initial 
distribution of objects was in near-circular, coplanar orbits in an annulus 
with a between 29 and 35 AU (designed to yield mostly Pluto-like orbits 
locked in the 3:2 orbital resonance with Neptune and with e between 0.2 and 
0.3). The time scale of orbital migration of the planets in this run was 
r=2X 10 6 yr. 

maintaining long-term orbital stability. The large inclinations 
and the argument-of-perihelion librations found in the work 
presented here augment the previous results and enhance the 
plausibility of the resonance capture theory for the origin of 
Pluto’s orbit. However, it is also clear that there is a rather 
small probability for obtaining inclinations as high as Plu- 
to’s. 

From the observed characteristics of Pluto’s orbit, I have 
shown that it is possible to infer the “initial” locations of 
both Neptune and Pluto and thus obtain a lower limit of 
about 5 AU for the magnitude of Neptune’s orbital expan- 
sion. (With further modeling, one can anticipate that it will 
be possible to also self-consistently infer the “initial” orbits 
of the other giant planets.) These inferences are of consider- 
able import for many aspects of planet formation. Some ex- 
amples are: the provenance of the outer planets may derive 
from a different and perhaps larger range of heliocentric dis- 
tance than the immediate vicinity of their current orbits; Ju- 
piter’s orbital migration could have influenced the planetesi- 
mal dynamics and accumulation process in the region of the 
asteroid belt; the capture of irregular satellites may have 
been significantly aided by the radial migration of the 
(proto-)giant planets. However, a discussion of these points 
is beyond the scope of the present work. 

In this paper, I have followed this “resonance capture” 


Fig. 6. The same as Fig. 5 but for Run 4 in which the time scale of orbital 
migration of the planets was t=4X10 6 yr. 

scenario for the origin of Pluto’s orbit to its logical next step 
in studying its implications for the architecture of the solar 
system beyond Neptune, i.e., the “Kuiper Belt” of comets 
approximately between Neptune’s orbit and 50 AU The nu- 
merical experiments reported here indicate that the dynami- 
cal structure of the Kuiper Belt is dominated by concentra- 
tions of objects trapped in orbital resonances with Neptune, 
particularly at the 3:2 and the 2:1 resonances. These resonant 
objects move on highly eccentric orbits, with a significant 
fraction on Neptune-crossing orbits; the inclinations of most 
of the objects remain low (less than 10°), but a small fraction 
(up to ~10%) are in the 15°-20° range. Libration of the 
argument of perihelion (about ±90°) is not an uncommon 
occurrence among the resonant objects. 

The numerical results also show that the inclination dis- 
tribution is sensitive to the rate of orbital evolution of the 
giant planets: longer time scale of the orbit evolution is cor- 
related with higher inclinations (cf. the third panel in Figs. 
5-7). The time scales used in the numerical experiments 
were chosen primarily for computational convenience, and 
ranged from 2 to 10 million years. As the libration period for 
the argument-of-perihelion libration (which is correlated 
with inclination excitation) is several million years, it is clear 
that the time scale of this resonance is comparable to the 
time scale of orbital evolution of the giant planets in these 
numerical experiments. Thus, the sensitivity of the inclina- 
tion distribution to the rate of orbital evolution of the planets 
is not surprising as the evolution is not “adiabatic” on this 
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Fig. 7. The same as Figs. 5 and 6, but for Run 7 in which the time scale of 
orbital migration of the planets was t= 1X10 7 yr. Note that there is little 
variation in the eccentricity distributions in Figs. 5-7, but the inclination 
distribution is progressively broader. 

time scale. It is worth noting that the numerical experiments 
of Fernandez & Ip (1984) indicate orbital evolution time 
scales of several tens of millions of years, and I expect that 
more realistic modeling in the future will help determine this 
time scale better. Comparison of these inclination distribu- 
tions with observations (when statistically significant num- 
bers of observations of trans-Neptunian objects and their or- 
bital elements are available) may provide a diagnostic of the 
rate of radial migration of the planets during the late stages 
of their evolution. 

Of course, the integrations reported here were only 
2X10 7 -10 8 yr long, and one might ask if the subsequent 
nearly conservative evolution over 10 9 yr time scales might 
not change the profile of the surviving objects. I have ex- 
tended the orbit integrations to 1 Byr in a few individual 
cases and have found that the resonance protection remained 
in place over this longer time period. Although suggestive, 
this is of course not proof of orbital stability over the 4.5 Byr 
age of the solar system. A slow leakage of objects out of the 
orbital resonances over this time scale — either due to colli- 
sional evolution, or possibly due to long-term purely gravi- 
tational effects, or both — cannot be ruled out, and may be 
necessary to supply the short-period comet population. In 
this context, it is also worthwhile to note that resonance cap- 
ture is not 100% efficient: a small fraction of the original 
trans-Neptunian population survives the “resonance sweep- 
ing” in its primordial nonresonant low-eccentricity, low- 


inclination orbits, and may also contribute to the fit x of 
short-period comets. A quantitative study of this poirit re- 
mains to be done. 

5.7 Comparison with Previous Theoretical Studies 
of the Kuiper Belt 

There have been two previous theoretical investigations 
of the dynamics of small bodies in the outer solar system. 
Holman & Wisdom (1993) and Levison & Duncan (1993) 
have studied test particle orbital stability on 10 9 yr time 
scales. For test particles in nearly circular and low- 
inclination orbits beyond Neptune, these studies found or- 
bital instability on short time scales (<10 7 yr) interior to 
33-34 AU, an intricate structure of interspersed regions of 
stability and instability in the semimajor axes range of 34 to 
43 AU, and substantially stable orbits beyond 43 AU The 
intricacy of the dynamical structure appears to be particularly 
acute near orbital resonance locations, but as these features 
are not analyzed in detail in either paper, it is difficult to 
draw clear conclusions about the relative number density of 
Kuiper Belt objects that might be expected in these regions 
at the present epoch. Holman and Wisdom also found sub- 
stantial “bumps” in the eccentricities and inclinations near 
the 3:2 resonance and in the eccentricity at the 2:1 resonance, 
indicating potential sources of planet-crossing cometary or- 
bits in these zones. Both studies indicate that once an object 
became Neptune crossing, it typically suffered a close en- 
counter with the planet in short order (although a few excep- 
tional cases of long-lived Neptune crossers were detected). 
In other words, Neptune-crossing orbits do not typically en- 
joy long-term stability due to resonance protection even 
when they originate close to resonance locations. This s not 
surprising, for the stable orbits near Neptune resonances ex- 
ist only in exceedingly narrow zones in phase space, ind a 
random sampling of initial conditions would have low prob- 
ability of hitting these zones. 

The resonance sweeping mechanism studied ir the 
present paper predicts the inner edge of stability at approxi- 
mately 36 AU in the semimajor axis, and clear and strong 
enhancements in the stability and relative population density 
of objects in orbital resonances with Neptune. In effect, the 
extensive nonresonant regions of space are depleted of ob- 
jects as those objects are swept into resonant orbits. In spite 
of the low phase space volume represented by the resonance 
libration regions, the dissipative nature of the dynamica evo- 
lution in this model very strongly enhances the occupancy 
rate of the resonance zones. This model predicts somewhat 
greater “dynamical erosion” in the Kuiper Belt than the 
above-mentioned studies of Holman and Wisdom and Levi- 
son and Duncan, particularly in the regions in between the 
resonances. The high eccentricities of the resonance-trapped 
objects make a large fraction of them Neptune crossirg yet 
dynamically long-lived. 

5.2 Implications For Observational Surveys 
of the Kuiper Belt 

The chief dynamical feature of the resonance-locked or- 
bits is that the longitude of perihelion avoids Neptune. (This 
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Table 1. The first column is magnitude, m\ s(peri) and s(aph) are the mini- 
mum radii (in km) of objects of magnitude less than m at perihelion and 
aphelion, respectively, in orbits of eccentricity 0.2. Columns 2 and 3 are for 
objects in the 3:2 Neptune resonance (semimajor axis =39.4 AU); columns 
4 and 5 are for objects in the 2:1 Neptune resonance (semimajor axis=47.8 
AU). We assume a mean geometric albedo A = 0.1. 


3:2 resonance 2:1 resonance 

m s (peri) s(aph) s(peri) s(aph) 


17 

690 

1570 

1020 

2324 

20 

173 

395 

256 

584 

22 

69 

157 

102 

232 

25 

17 

39 

26 

58 


fact can be inferred from the librations of the resonance 
angle.) For example, in the case of the 3:2 resonance, the 
resonance angle librates about 180°, and its maximum libra- 
tion amplitude (for stable orbits) is about 90°; the longitude 
of perihelion librates about a center ±90° away from Nep- 
tune’s mean longitude. The situation is more complicated for 
the 2:1 resonance, as the location of the center of libration is 
a strong function of the orbital eccentricity; the center of the 
perihelion oscillation (relative to Neptune’s mean longitude) 
is at 180° for very small eccentricities; for larger eccentrici- 
ties (exceeding —0.03), the libration center bifurcates and 
the two new libration centers drift away (symmetrically) 
from 180° to nearly +45° and —45° for eccentricities near 
0.3. These properties of the resonant orbits will be discussed 
in further detail in a future publication. 

If most Kuiper Belt objects do indeed exist primarily in 
two resonance bands (the 3:2 and the 2:1), are there particu- 
lar observational strategies that might improve the discovery 
rate of these objects? As the results in Sec. 4 indicate, the 
inclinations of these objects do not exceed 20° (more than 
80% of them remain less than 10°). Therefore, sky surveys 
within ±20° of the ecliptic are likely to yield most of the 
detectable objects. 2 In a magnitude-limited sky survey, the 
minimum size of a detectable outer solar system object will 
vary with its distance and its albedo. Assuming that observa- 
tions are made at opposition (geocentric distance 
= heliocentric distance -1 AU), the size-magnitude- 
heliocentric distance relationship is given by 

s~ 1803(A/0.1 )“ 1/2 r(r — 1 )10 _0 2m , (8) 

where s is the object radius in km, A is the geometric albedo, 
r is the heliocentric distance in AU, and m is the magnitude. 
For a typical orbital eccentricity of —0.2, the perihelion to 
aphelion distance varies from 31 to 47 AU in the 3:2 reso- 
nance, and from 38 to 58 AU in the 2:1 resonance. The 
minimum radii of detectable objects at perihelion and ap- 
helion in these resonant orbits are listed in Table 1 for several 
limiting magnitudes. For a given magnitude and albedo, the 
minimum radius of a detectable object varies by a factor of 


2 It should be noted that about a third of the resonant trapped objects also 
exhibit argument-of-perihelion libration; this means that at perihelion (when 
they are brightest and thus most likely to be detected), such objects would 
also be near their greatest ecliptic latitude (either above or below the eclip- 
tic); observational searches that are too narrowly confined to the ecliptic 
may not detect these objects. 


approximately — (1 T e) 2 /(l - e) 2 from perihelion to ap- 
helion. For e = 0.2, this factor is —2.3. For illustrative pur- 
poses, let us assume that the Kuiper Belt objects have a 
power law size distribution, n(s)ds^s~ q ds , with index 
q = 2 (cf. Jewitt & Luu 1995). Then, the number of 
resonance-trapped objects that are potentially detectable at 
perihelion is a factor —2.3 greater than those detectable at 
aphelion. Recall that the perihelia of the 3:2 resonant objects 
librate near ±90° from Neptune. The 2:1 resonant objects 
(with e — 0.2) have perihelia librating near +60° or —60° 
from Neptune. Assuming that, at any given epoch, the reso- 
nant objects are uniformly distributed in ecliptic longitude, 
those objects at longitudes defined by the two quadrants 
(45°, 135°) and ( — 135°,— 45°) relative to Neptune would be 
at or close to perihelion. Therefore, without going through a 
detailed statistical analysis (to take account of factors such as 
the distributions in eccentricity and resonance libration am- 
plitude), a crude estimate is that discovery rates in these 
regions of the sky would be greater by a factor of about 2 
compared to the two quadrants outside these ranges of lon- 
gitude. With Neptune currently near 290° ecliptic longitude, 
the favored quadrants correspond to ecliptic longitude ranges 
(335°, 65°) and (155°, 245°). 

It is interesting to consider the detectability of a Pluto-like 
object in the 2:1 resonance with Neptune. Such an object 
(with radius 5—1000 km, mean geometric albedo A— 0.5, 
and orbital eccentricity 0.25) would vary in brightness from 
m— 15.0 tom— 17.3 from perihelion to aphelion. Of all the 
sky surveys for outer solar system objects to date, Tom- 
baugh’s (1961) search is the most likely to have detected 
such a body. This survey covered all longitudes within ap- 
proximately ±15° of the ecliptic, with a limiting magnitude 
of about 16.5. Thus, a Pluto-like object in the 2:1 resonance 
would have been detected if it were at a heliocentric distance 
less than 49 AU In a 2:1 resonant orbit (semimajor axis 
=47.8 AU) with e = 0.25, an object spends only about 35% 
of its orbital period at heliocentric distances less than 49 AU 
Therefore, in Tombaugh’s survey, there was a —35% chance 
of detecting such an object. The more recent surveys for 
outer solar system objects (Kowal 1989; Luu & Jewitt 1988; 
Levison & Duncan 1990; Jewitt & Luu 1995) all had limit- 
ing magnitudes exceeding 17.3, but they also had much 
smaller sky coverage; detection probability is much smaller 
in these surveys for the latter reason. Therefore, a Pluto-like 
object in the 2:1 resonance cannot yet be ruled out. 

Observational surveys of the outer planetary system have 
recently reported the detection of several objects —100 km in 
size which are possibly the larger members of the Kuiper 
Belt (Jewitt & Luu 1995), and it appears likely that even 
greater numbers of detections will be forthcoming in the near 
future. As the present work was in progress, Marsden (1994a, 
b) has reported on the possibility that several of the newly 
discovered objects may be in Pluto-like orbits, locked in the 
3:2 resonance with Neptune. If this is confirmed, it would 
provide further corroboration for the ’‘resonance sweeping” 
scenario. It is my hope that the present paper will contribute 
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to the acceleration of these observational detections and their 
interpretation within models of the formation and evolution 
of the solar system. 

I am grateful to Scott Tremaine for a critical reading of 
this paper. This research was done while the author was a 


Staff Scientist at the Lunar and Planetary Institute which is 
operated by the Universities Space Research Association un- 
der Contract No. NASW-4574 with the National Aeronautics 
and Space Administration. This paper is Lunar and Planetary 
Institute Contribution No. 860. 
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